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Abstract 

The aim of these lectures is to show that the methods of classical Hamil- 
tonian mechanics can be profitably used to solve certain classes of non- 
linear partial differential equations. The prototype of these equations is 
the well-known Korteweg-de Vries (KdV) equation. 
In these lectures we touch the following subjects: 

i) the birth and the role of the method of Poisson pairs inside the 
theory of the KdV equation; 

the theoretical basis of the method of Poisson pairs; 
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the Gel'fand-Zakharevich theory of integrable systems on bi-Ham- 
iltonian manifolds; 

the Hamiltonian interpretation of the Sato picture of the KdV flows 
and of its linearization on an infinite-dimensional Grassmannian 
manifold. 

the reduction technique(s) and its use to construct classes of solu- 
tions; 

the role of the technique of separation of variables in the study of 
the reduced systems; 

some relations intertwining the method of Poisson pairs with the 
method of Lax pairs. 



x Work partially supported by the M.U.R.S.T. and the G.N.F.M. of the Italian C.N.R. 
Lectures given by F. Magri at the 1999 CIME course "Direct and Inverse Methods in Solving 
Nonlinear Evolution Equations " Cetraro, (Italy) September 1999. 
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1 Introduction: The tensorial approach and the 
birth of the method of Poisson pairs 

This lecture is an introduction to the Hamiltonian analysis of PDEs form an "ex- 
perimental" point of view. This means that we are more concerned in unveiling 
the spirit of the method than in working out the theoretical details. Therefore 
the style of the exposition will be informal, and proofs will be mainly omitted. 
We shall follow, step by step, the birth and the evolution of the Hamiltonian 
analysis of the KdV equation 

u t = ^{u xxx - Quu x ) , (1.1) 
from its "infancy" to the final representation of the KdV flow as a linear flow 



on an infinite-dimensional Grassmannian due to Sato ||27|| . The route is long 
and demanding. Therefore the exposition is divided in two parts, to be carried 
out in this and in the fourth lecture (see Section |]). Here our primary aim is 
to show the birth of the method of Poisson pairs. It is reached by means of a 
suitable use of the well-known methods of tensor analysis. We proceed in three 
steps. First, by using the transformation laws of vector fields, we construct 
the Miura map and the so called modified KdV equation (mKdV). This result 
leads quite simply to the theory of (elementary) Darboux transformations and 
to the concept of Poisson pair. Indeed, a peculiarity of mKdV is to possess 
an elementary Hamiltonian structure. By means of the transformation law 
of Poisson bivectors, we are then able to transplant this structure to the KdV 
equation, unraveling its "bi-Hamiltonian structure" . This structure can be used 
in turn to define the concept of Lenard chain and to plunge the KdV equation 
into the "KdV hierarchy" . This step is rather important from the point of view 
of finding classes of solutions to the KdV equation. Indeed the hierarchy is a 
powerful instrument to construct finite-dimensional invariant submanifolds of 
the equation and, therefore, finite-dimensional reductions of the KdV equation. 
The study of this process of restriction and of its use to construct solutions will 
be one of the two leading themes of these lectures. It is intimately related to the 
theory of separation of variables dealt with in the last two lectures. The second 
theme is that of the linearization of the full KdV flow on the infinite-dimensional 
Sato Grassmannian. The starting point of this process is surprisingly simple, 
and once again based on a simple procedure of tensor calculus. By means of the 
transformation laws of one-forms, we pull back the KdV hierarchy from its phase 
space onto the phase space of the mKdV equation. In this way we obtain the 
"mKdV hierarchy" . In the fourth lecture we shall show that this hierarchy can 
be written as a flow on an infinite-dimensional Grassmann manifold, and that 
this flow can be linearized by means of a (generalized) Darboux transformation. 

1.1 The Miura map and the KdV equation 

As an effective way of probing the properties of equation (|1 . 1|) we follow the 
tensorial approach. Accordingly, we regard equation ( [Lip as the definition of a 
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vector field 



Ut — X(u,u x ,u xx ,u xxx ) (1-2) 

on a suitable function space, and we investigate how it transforms under a point 
transformation in this space. Since our "coordinate" u is a function and not 
simply a number, we are allowed to consider transformations of coordinates 
depending also on the derivatives of the new coordinate function of the typeQ 

u = <S>(h,h x ). (1.3) 

We ask whether there exists a transformed vector field 

h t = Y(h,h x ,h xx ,h xxx ) (1.4) 

related to the KdV equation according to the transformation law for vector 
fields, 

X($(h)) = & h -Y(h), (1.5) 

where Q' h is the (Frechet) derivative of the operator $ defining the transfor- 
mation. This condition gives rise to a (generally speaking) over-determined 
system of partial differential equations on the unknown functions $(/i, h x ) and 
Y(h, h x , h xx , h xxx ). In the specific example the over-determined system can be 
solved. Apart form the trivial solution u = h x , we find the Miura transformation 



u = h x + h 2 -X, (1.6) 

depending on an arbitrary parameter A. The transformed equation is the mod- 
ified KdV equation: 

h t = ~{h xxx - 6h 2 h x + 6Xh x ) . (1.7) 



Exercise 1.1 Work out in detail the transformation law ( |1.5|) , checking that 
X, $, and Y, defined respectively by equations (|1.1| ), (|1.6|) and ( |1.7| ) do satisfy 
equation (|1.5|) . □ 

The above result is plenty of consequences. The first one is a simple method 
for constructing solutions of the KdV equation. It is called the method of 
(elementary) Darboux transformations p2| . It rests on the remark that the 
mKdV equation ( |L~T| ) admits the discrete symmetry 

h^h' = -h. (1.8) 

Let us exploit this property to construct the well-known one-soliton solution of 
the KdV equation. We notice that the point u = is a very simple (singular) 

2 For further details on these kind of transformations, see M. 
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invariant submanifold of the KdV equation. Its inverse image under the Miura 
transformation is the 1 -dimensional submanifold S\ formed by the solutions of 
the special Riccati equation 

h x + h 2 = X. (1.9) 

This submanifold, in its turn, is invariant with respect to equation ( |1.7|) . A 
straightforward computation shows that, on this submanifold, 

\{h xxx - 6h 2 h x + 6\h x ) = Xh x . (1.10) 

Therefore, on Si the mKdV equation takes the simple form h t = Xh x . Solving 
the first order system formed by this equation and the the Riccati equation 
h x + h 2 = A, and setting A = z 2 , we find the general solution 

h(x,t) = z tanh(zx + zH + c) (1.11) 

of the mKdV equation on the invariant submanifold S±. At this point we use 
the symmetry property and the Miura map. By the symmetry property (|1.8| ) 
the function —h(x,t) is a new solution of the modified equation, and by the 
Miura map the function 

u '(x, t) = -h x + h 2 - z 2 = 2z 2 sech 2 (zx + z 3 t + c) (1.12) 

is a new solution of the KdV equation. It is called the one soliton solution^. It 
can also be interpreted in terms of invariant submanifolds. To this end, we have 
to notice that the Miura map ( |1.6| ) transforms the invariant submanifold Si of 



the modified equation into the submanifold formed by the solutions of the first 
order differential equation 

lf-l u l + u ^j+Xu 2 = 0. (1.13) 

As one can easily check, this set is preserved by the KdV equation, and therefore 
it is an invariant one-dimensional submanifold of the KdV equation, built up 
from the singular manifold u — 0. On this submanifold, the KdV equation 
takes the simple form ut = Xu x , and the flow can be integrated to recover the 

solution (|rrj). 

This example clearly shows that the Darboux transformations are a mech- 
anism to build invariant submanifolds of the KdV equation. Some of these 
submanifolds will be examined in great detail in the present lectures. The pur- 
pose is to show that the reduced equations on these submanifolds are classical 
Hamiltonian vector fields whose associated Hamilton-Jacobi equations can be 
solved by separation of variables. In this way, we hope, the interest of the 
Hamiltonian analysis of the KdV equations can better emerge. 



3 For a very nice account of the origin and of the properties of the KdV equation and of 
other soliton equations and their solutions, see, e.g., p3. 
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1.2 Poisson pairs and the KdV hierarchy 

We shall now examine a more deep and far reaching property of the Miura map. 
It is connected with the concept of Hamiltonian vector field. From Analytical 
Mechanics, we know that the Hamiltonian vector fields are the images of exact 
one-forms through a suitable linear map, associated with a so-called Poisson 
bivector. We shall formally define these notions in the next lecture. These defi- 
nition can be easily extended to vector fields on infinite-dimensional manifolds. 
Let us give an example, by showing that the mKdV equation is a Hamiltonian 
vector field. This requires a series of three consecutive remarks. First we notice 
that equation (|1.7Q can be factorized as 



1 3 

-h T x h 3 + 3Xh 

2 2 



;i.i4) 



Then, we notice that the linear operator in the first bracket, \d X) is a constant 
skewsymmetric operator which we can recognize as a Poisson bivector. Finally, 
we notice that in the differential polynomial appearing in the second bracket in 
the right hand side of equation ( |1.14|) , we can easily recognize an exact one- 
form. Indeed, 

-h xx - -h 3 + 3Xh] hdx = ^- [ (--hi - -h 4 + -Xh 2 ] dx (1.15) 
2 2 J dt J \ 2 8 2 y 

for any tangent vector h. These statements are true under suitable boundary 
conditions, as explained in, e.g., ||. Here and in the rest of these lectures we 
will tacitly use periodic boundary conditions. 

The Hamiltonian character of the mKdV equation is obviously independent 
of the existence of the Miura map. However, this map finely combines this 
property from the point of view of tensor analysis. Let us recall that a Pois- 
son bivector is a skewsymmetric linear map from the cotangent to the tangent 
spaces satisfying a suitable differential condition (see Lecture |2|). It obeys the 
transformation law 

Q m = <!>' h P h & h * (1.16) 

under a change of coordinates (or a map between two different manifolds). In 
this formula the point transformation is denoted (in operator form) by u = 
$(/i). The operators <&' h and are the Frechet derivative of $ and its adjoint 
operator. The symbols Ph and Q u denote the Poisson bivectors in the space of 
the functions h and u, respectively. Since the Miura map u = h x + h 2 — X is not 
invertible, it is rather nontrivial that there exists a Poisson bivector Q u , on the 
phase space of the KdV equation, which is ^-related (in the sense of equation 
(|1.16|) ) to the Poisson bivector P h = jd x associated with the modified equation. 
Surprisingly, this is the case. One can check that the operator Q u is defined by 

Q u = ~\d xxx + 2(u + X)d x + u x . (1.17) 
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Exercise 1.2 Verify the above claim by computing the product (in the appro- 
priate order) of the operators &' h = d x + h, P h = \d x , and = — d x + h, and 
by expressing the results in term of u = h x + h 2 — A. □ 

This exercise shows that the Miura map is a peculiar Poisson map. Since it 
depends on the parameter A, the final result is that the phase space of the KdV 
equation is endowed with a one-parameter family of Poisson bivectors, 

Qa = Qi-AQ , (1.18) 

which we call a Poisson pencil. The operators (Qi,Qo) defining the pencil are 
said to form a Poisson pair, a concept to be systematically investigated in the 
next lecture. 

These operators enjoy a number of interesting properties, and define new 
geometrical structures associated with the equation. One of the simplest but 
far-reaching is the concept of Lenard chain. The idea is to use the pair of 
bivectors to define a recursion relation on one-forms: 



QoQj+i = Qiotj • (1-19) 

In the applications a certain care must be taken in dealing with this recursion 
relation, since it does not define uniquely the forms otj (the operator Q is 
seldom invertible). Furthermore, it is still less apparent that it can be solved 
in the class of exact one-forms. However, in the KdV case we bonafide proceed 
and we find 

a Q — 1 

1 

«2 = -{-u xx + 3u ) 

"3 = ^(~ 10m3 + 10 uu xx - u xxxx + 5u 2 x ) 

as first terms of the recurrence. The next step is to consider the associated vector 
fields (the meaning of numbering them with odd integers will be explained in 
Lecture 4): 

du n n 

— = Qia = Qotti = u x 

at i 

du ! 

— = Q iai = Q a 2 = -( Uxxx - Q UUx ) (1.21) 
du i 

— = Qxa 2 = Qoa 3 = —(u xxxxx - I0uu xxx - 20u x u xx + 30u 2 u x ) . 
ot 5 16 

They are the first members of the KdV hierarchy. In the fourth lecture, we shall 
show that it is a special instance of a general concept, the Gel'fand-Zakharevich 
hierarchy associated with any Poisson pencil of a suitable class. 
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1.3 Invariant submanifolds and reduced equations 

The introduction of the KdV hierarchy has important consequences on the prob- 
lem of constructing solutions of the KdV equation. The hierarchy is indeed a 
basic supply of invariant submanifolds of the KdV equation. This is due to the 
property of the vector fields of the hierarchy to commute among themselves. 
From this property it follows that the set of singular point of any linear com- 
bination (with constant coefficients) of the vector fields of the hierarchy is a 
finite-dimensional invariant submanifold of the KdV flow. These submanifold 
can be usefully exploited to construct classes of solutions of this equation. 

As a first example of this technique we consider the submanifold defined by 
the condition 

that is, the submanifold where the second vector field of the hierarchy is a 
constant multiple of the first one. It is formed by the solutions of the third 
order differential equation 

j(u xxx - 6uu x ) - Xu x = . (1.23) 

Therefore it is a three dimensional manifold, which we denote by M3. We can 
use as coordinates on M 3 the values of the function u and its derivatives u x and 
u xx at any point xq. To avoid cumbersome notations, we will continue to denote 
these three numbers with the same symbols, u,u x ,u xx , but the reader should 
be aware of this subtlety. To perform the reduction of the first equation of the 
hierarchy ( 1.21|) on M 3 , we consider the first two differential consequences of 



the equation = u x and we use the constraint ( 1.22j) to eliminate the third 
order derivative. We obtain the system 

du du x du xx 

77T = Ux, -77- = u xx , — — = 6uu x + 4Xu x . (1.24) 
ati oti oti 

We call X\ the vector field defined by these equations on M3. It shares many 
of the properties of the KdV equation. For instance, it is related to a Poisson 
pair. The simplest way to display this property is to remark that X\ possesses 
two integrals of motion, 

Hq = u xx — 3u 2 — AXu 

1 2 , 2 (1-25) 

H x = —ul + u 3 + 2Xu 2 + uH . 

Then we notice that on M3 there exists a unique Poisson bracket {•, -} with 
the following two properties: 

i) The function Hq is a Casimir function, that is, {F, Hq} = for every 
smooth function F on M3. 
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ii) X\ is the Hamiltonian vector field associated with the function Hi. 

Such a Poisson bracket {•, -} is defined by the relations 

{u,u x } Q = -l, {u,u xx } = 0, {u x ,u xx } = 6u + 4A . (1.26) 

Similarly, one can notice that on M3 there exists a unique Poisson bracket {•, -}i 
with the following "dual" properties: 

i') The function Hi is a Casimir function, that is, {F,Hi} 1 = for every 
smooth function F on M 3 . 

ii') X\ is the Hamiltonian vector field associated with the function Hq. 

This second Poisson bracket {•, •} 1 is defined by the relations 

{u,u x } 1 = u, {u,u xx } 1 = u x , {u x , u xx } 1 = u xx — u(6u + 4A) . (1.27) 



Exercise 1.3 Verify the stated properties of the Poisson pair (|1.26|) and ( |1.27|) . 



□ 

We now exploit the previous remarks to understand the geometry of the flow 
associated with X\. First we use the Hamiltonian representation 

dF 

— = X 1 (F) = {F,H 1 } . (1.28) 

It entails that the level surfaces of the Casimir function H are two-dimensional 
(symplectic) manifolds to which Xi is tangent. Let us pick up any of these 
symplectic leaves, for instance the one passing through the origin u = 0, u x = 
0, u xx = 0. Let us call S 2 this leaf. The coordinates (u, u x ) are canonical co- 
ordinates on S 2 - The level curves of the Hamiltonian Hi define a Lagrangian 
foliation of S2. Our problem is to find the flow of the vector field X\ along 
these Lagrangian submanifold. We have already given the solution of this prob- 
lem in the particular case of the Lagrangian submanifold passing through the 
origin. This submanifold is the one-dimensional invariant submanifold ( 1.13Q 



previously discussed in connection with Darboux transformations. The relative 
flow is the one-soliton solution to KdV. To deal with the generic Lagrangian 
submanifolds on an equal footing, it is useful to change strategy and to use the 
Hamilton-Jacobi equation 

H 1 (u,^-) = E. (1.29) 
du 

In this rather simple example, there is almost nothing to say about this equation 
(which is obviously solvable), and the second Poisson bracket ( |1.26| ) seems not 
to play any role in the theory. 

This (wrong!) impression is promptly corrected by the study of a more 
elaborated example. Let us consider the five-dimensional submanifold M5 of 
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the singular points of the third vector field of the KdV hierarchy. It is defined 
by the equation 



u, 



10uu xxx — 20u x u xx + 30u 2 u x = . (1.30) 



On this submanifold we can consider the restrictions of the first two vector fields 
of the same hierarchy. To compute the reduced equation we proceed as before. 
We regard the Cauchy data (u,u x ,u xx ,u xxx ,u xxxx ) as coordinates on M 5 . then 

i • i • • dtJj QUi^fc ^^-^XX ^'^■'XXX ^'^■'XXXX i n • i 

we compute the time derivatives — — , — — , — — , — — , — by taking the 

ot\ at\ oti oti oti 

differential consequences of = u x , and by using ( |1.30|) and its differential 

consequences as a constraint to eliminate all the derivatives of degree higher 
than four. We obtain the equations 



du 
df x 
du x 

du xx 
~dtT 

du xxx 

dt 1 

d^Tcr.cr.cr. 



U.7 



(1.31) 



10uu xxx + 20u x u xx — 30u 2 u x 



In the same way, for the reduction of the KdV equation, we get 



du 1 



- m ttt — 6mm 



dt 3 4 
du T 1 

XXXX U U U'u 



^XXXX SUUtt 6li™ 



dt 3 4 

du rr 1 



{Auu xxxx + 2u x u xx — 30u 2 u x ) 



dt 3 4 " ' (1-32) 

r^ xx = ^(2«L + Qu x u xxx + 4uu xxxx - 60uu 2 x - 30u 2 u xx ) 
du 1 

— ^ xxx = ~(WuxxU xxx + Wu x u xxxx — 120n 3 n x — 100uu x u xx 
dt 3 4 

+ 10u 2 u xxx - 60ul) . 



Exercise 1.4 Verify the previous computations. □ 

To find the corresponding solutions of the KdV equation, regarded as a partial 
differential equation in x and t, we have to: 



1. Construct a common solution to the ordinary differential equations ( 1.31]) 
and ([TBI ); 
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2. Consider the first component u(ti,ts) of such a solution; 



3. Set t\ = x and £3 = t. 

The function u(x,t) obtained in this way is the solution we were looking for. 
What makes this procedure worth of interest is that the ODEs ( 1.31 ) — ( pL.32| ) 



can be solved by means of a variety of methods. In particular, they can be 
solved by means of the method of separation of variables^. It can be shown 
that they are rather special equations: They are Hamiltonian with respect to 
a Poisson pair; this Poisson pair allows to foliate the manifold M5 into four- 
dimensional symplectic leaves with special properties; each symplectic leaf 64 
carries a Lagrangian foliation to which the vector fields (|1.31|) and (|1.32|) are 



tangent; the Poisson pair defines a special set of coordinates on each S4; in these 
coordinates the Hamilton-Jacobi equations associated with the Hamiltonian 
equations ( |1.31| ) and ( |1.32| ) can be simultaneously solved by additive separation 



of variables. Most of these properties will be proved in the next lecture. 
1.4 The modified KdV hierarchy 

We leave for the moment the theme of the reduction, and come back to the 
KdV hierarchy in its general form. We notice that the first equations (|1.21| ) can 
also be written in the form 



du 
dh 
du 

-wr = {Qi- \Qo){\ 2 at Q + A«i + a 2 ) 
ot 5 



(Qi - \Qo)a Q 

(Qi-XQ )(Xa + ai) (1.33) 



This representation shows that these equations are Hamiltonian with respect to 
the whole Poisson pencil. This elementary property can be exploited to simply 
construct the modified KdV hierarchy. Let us write in general 

°" ={Q 1 + \Q )aW(\), (1.34) 



dt2j+i 
where 

«W(A) = \ j a + A j_1 ai + • • • + a d . (1.35) 



By means of the Miura map (|1.6f) we pull-back the one-forms a^' to one- 
form (3^ defined on the phase space of the modified equation according to the 
transformation law of one-forms, 

(3(h) = & h *a(Hh)) ■ (1-36) 



4 The fact that the stationary reductions of KdV can be solved by separation of variables 
is well-known (see, e.g., |)|). This classical method has recently found a lot of interesting new 
applications, as shown in the survey pq|. 
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We then define the equations 

dh 



2j+l 



PkP U) (\) ■ (1-37) 



They are ^-related to the corresponding equations of the KdV hierarchy, exactly 
as the mKdV equation is ^-related to the KdV equation. Indeed, 

0(1 = Wa^^W; A ) = Qu^\u; A) . (1.38) 



It is therefore natural to call equations (11.371 ) the modified KdV hierarchy. By 
using the explicit form of the operators Ph and $^*, it is easy to check that the 
modified hierarchy is defined by the conservation laws 

<>k - = d x H& +1 \ (1.39) 



dt2j+i 



where 



H^ +1) = —a^ + a^h . (1.40) 



Exercise 1.5 Write down explicitly the first three equations of the modified 
hierarchy. □ 

The above formulas are basic in the Sato approach. In the fourth lecture, after 
a more accurate analysis of the Hamiltonian structure of the KdV equations, we 
shall be led to consider the currents H^ +1 ' as defining a point of an infinite- 
dimensional Grassmannian. This point evolves in time as the point u moves 
according to the KdV equation. We shall determine the equation of motion of 
the currents H^ 2 ^ +l \ They define a "bigger" hierarchy called the Central System. 
This system contains the KdV hierarchy as a particular reduction. It enjoys the 
property of being linearizable. In this way, by a continuous process of extension 
motivated by the Hamiltonian structure of the equations (from the single KdV 
equation to the KdV hierarchy and to the Central System), we arrive at the 
result that the KdV flow can be linearized. At this point the following picture 
of the possible strategies for solving the KdV equations emerges: either we pass 
to the Sato infinite-dimensional Grassmannian and we use a linearization tech- 
nique, or we restrict the equation to a finite-dimensional invariant submanifold 
and we use a technique of separation of variables. The two strategies comple- 
ment themselves rather well. Our attitude is to see the Grassmannian picture 
as a compact way of defining the equations, and the "reductionist" picture as 
an effective way for finding interesting classes of solutions. 
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The plan of the lectures 



This is the web of ideas which we would like to make more precise in the following 
lectures. As cornerstone of our presentation we choose the concept of Poisson 
pairs. In the second lecture, we develop the theory of these pairs up to the point 
of giving a sound basis to the concept of Lenard chain. In the third lecture we 
exhibit a first class of examples, and we explain a reduction technique allowing 
to construct the Poisson pairs of the reduced flows. In the fourth lecture we give 
a second look at the KdV theory, and we explain the reasons which, according 
to the Hamiltonian standpoint, suggest to pass on the infinite-dimensional Sato 
Grassmannian. In the fifth lecture we better explore the relation between the 
two strategies, and we touch the concept of Lax representation. Finally, the last 
lecture is devoted to the method of separation of variables. The purpose is to 
show how the geometry of the reduced Poisson pairs can be used to define the 
separation coordinates. 
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2 The method of Poisson pairs 



In the previous lecture we have outlined the birth of the method of Poisson 
pairs and its main purpose: To define integrable hierarchies of vector fields. In 
this lecture we dwell on the theoretical basis of this construction presenting the 
concept of Gel'fand-Zakharevich system. 

The starting point is the notion of Poisson manifold. A manifold is said 
to be a Poisson manifolof] if a composition law on scalar functions has been 
defined obeying the usual properties of a Poisson bracket: bilinearity, skewsym- 
metry, Jacobi identity and Leibnitz rule. The last condition means that that 
the Poisson bracket is a derivation in each entry: 

{fg,h} = {f,h}g + f{g,h} . (2.1) 

Therefore, by fixing the argument of one of the two entries and keeping free the 
remaining one, we obtain a vector field, 

X h = {;h}. (2.2) 

It is called the Hamiltonian vector field associated with the function h with 
respect to the given Poisson bracket. Due to the remaining conditions on the 
Poisson bracket, these vector fields are closed with respect to the commutator. 
They form a Lie algebra homeomorphic to the algebra of functions defined by 
the Poisson bracket: 

\ X Ji X g\ = X {f,g} ■ ( 2 - 3 ) 

Therefore a Poisson bracket on a manifold has a twofold role: it defines a Lie 
algebra structure on the ring of C°°-functions, and provides a representation of 
this algebra on the manifold by means of the Hamiltonian vector fields. 

Instead of working with the Poisson bracket, it is often suitable to work 
(especially in the infinite-dimensional case) with the associated Poisson tensor. 
It is the bivector field P on M defined by 

{f,g} = (df, Pdg) . (2.4) 

In local coordinates, its components P 3 (x 1 , . . . , x n ) are the Poisson brackets of 
the coordinate functions, 

P jk {x\... ,x n ) = {x 3 ,x k } . (2.5) 

By looking at this bivector field as a linear skewsymmetric map P : T*M — > 
TM, we can define the Hamiltonian vector fields Xf as the images through P 
of the exact one-forms, 

X f = Pdf. (2.6) 

In local coordinates this means 

*1 = Pik ^ ■ ( 2 - 7 ) 



5 The books (T§ and |T) are very good references for this topic. 
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Exercise 2.1 Show that the components of the Poisson tensor satisfy the cyclic 
condition 

p "-ir + pt '%r^%-)=o. (2.8) 

□ 



Exercise 2.2 Suppose that M is an affine space A. Call V the vector space 
associated with A. Define a bivector field on A as a mapping P : A x V* — > V 
which satisfies the skewsymmetry condition 

(a, P u (3) = -(J3, P u a) (2.9) 

for every pair of covector (a,/3) in V* and at each point u G A. Denote the 
directional derivative at the point u of the mapping u \— > P u a along the vector 
v by 

P^a;v) = j t P u+tv a\ t=0 . (2.10) 

Show that the bivector P is a Poisson bivector if and only if it satisfies the cyclic 
condition 

(a, P'M P ul )) + (/3, J*( 7 ; P u a)) + ( 7 , P>; P U P)) = • (2.11) 

□ 



Exercise 2.3 Check that the bivector Q\ of equation (|1.17|) , associated with 
the KdV equation, fulfills the conditions (p.9|) and (|2.11| ). □ 



No condition is usually imposed on the rank of the Poisson bracket, that is, 
on the dimension of the vector space spanned by the Hamiltonian vector fields 
at each point of the manifold. If these vector fields span the whole tangent 
space the bracket is said to be regular, and the manifold M turns out to be 
a symplectic manifold. Indeed there exists, in this case, a unique symplectic 
2-form oj such that 



{f,g} = u(X f ,X 9 ) . (2.12) 

More interesting is the case where the bracket is singular. In this case, the 
Hamiltonian vector fields span a proper distribution D on M. It is involutive 
but, generically, not of constant rank. Nonetheless, this distribution is com- 
pletely integrable: at each point there exists an integral submanifold of maximal 
dimension which is tangent to the distribution. These submanifolds are sym- 
plectic manifolds, and are called the symplectic leaves of the Poisson structure. 
The symplectic form is still defined by equation ( |2.12| ). Indeed, even if there 
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is not a 1-1 correspondence between (differentials of) functions and Hamilto- 
nian vector fields, this formula keeps its meaning, since the value of the Poisson 
bracket does not depend on the particular choice of the Hamiltonian function 
associated with a given Hamiltonian vector field. We arrive thus at the following 
conclusion: a Poisson manifold is either a symplectic manifold or a stratification 
of symplectic manifolds possibly of different dimensions. It can be proven that, 
in a sufficiently small open set where the rank of the Poisson tensor is con- 
stant, these symplectic manifolds are the level sets of some smooth functions 
Fx, . . . , Fk, whose differentials span the kernel of the Poisson tensor. They are 
called Casimir functions of P (see below). 

Exercise 2.4 Let {x\,x 2 , x 3 } be Cartesian coordinates in M = M 3 . Prove that 
the assignment 

{xi,x 2 } = x 3 , {x 1 ,x 3 } = -x 2 , {x 2 ,x 3 } = x 1 (2.13) 

defines a Poisson tensor on M. Find its Casimir function, and describe the 
symplectic foliation associated with it. □ 

After these brief preliminaries on Poisson manifolds as natural settings for the 
theory of Hamiltonian vector fields, we pass to the theory of bi-Hamiltonian 
manifolds. Our purpose is to provide evidence that they are a suitable setting 
for the theory of integrable Hamiltonian vector fields. The simplest connection 
between the theory of integrable Hamiltonian vector fields and the theory of bi- 



Hamiltonian manifold is given by the Gel'fand-Zakharevich (GZ) theorem JT3 
14|1 we shall discuss in this lecture. 



A bi-Hamiltonian manifold is a Poisson manifold endowed with a pair of 
compatible Poisson brackets. We shall denote these brackets with {f,g}o and 
{f,g}i- They are compatible if the Poisson pencil 

{f,g}x:={f,g}i-\{f,g}o (2.14) 

verifies the Jacobi identity for any value of the continuous (say, real) parameter 
A. By means of this concept we catch the main features of the situation first 
met in the KdV example of Lecture |I[ The new feature deserving attention is 
the dependence of the Poisson bracket ( [2.14 ) on the parameter A. It influences 



all the objects so far introduced on a Poisson manifold: Hamiltonian fields and 
symplectic foliation. In particular, this foliation changes with A. The useful 
idea is to extract from this moving foliation the invariant part. It is defined 
as the intersection of the symplectic leaves of the pencil when A ranges over 
ffi. U {oo}. The GZ theorem describes the structure of these intersections in 
particular cases. 

Let us suppose that the dimension of M is odd, dimM = 2n+ 1, and that 
the rank of the Poisson pencil is maximal. This means that the dimension of 
the characteristic distribution spanned by the Hamiltonian vector field is 2n for 
almost all the values of the parameter A, and almost everywhere on the manifold 
M. In this situation the generic symplectic leaf of the pencil has accordingly 
dimension 2n and the intersection of all the symplectic leaves are submanifolds 
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of dimension n. For brevity, we shall call this intersection the support of the 
pencil. The GZ theorem displays an important property of the leaves of the 
support of the pencil. 

Theorem 2.5 On a (2n+l) -dimensional bi-Hamiltonian manifold, whose Pois- 
son pencil has maximal rank, the leaves of the support are generically Lagrangian 
submanifolds of dimension n contained on each symplectic leaf of dimension 2n. 

This theorem contains two different statements. First of all it states that the 
dimension of the support is exactly half of the dimension of the generic sym- 
plectic leaf. It is the "hard" part of the theorem. Then it claims that the leaves 
of the support are Lagrangian submanifolds. Contrary to the appearances, this 
is the easiest part of the theorem, as we shall see. To better understand the 
content of the GZ theorem, we deem suitable to look at it from a different and, 
so to say, more constructive, point of view. It requires the use of the concept 
of Casimir function, defined as a function which commutes with all the other 
functions with respect to the Poisson bracket. Equivalently, it can be defined 
as a function whose differential belongs to the kernel of the Poisson tensor, i.e., 
a function generating the null vector field. In the case of a Poisson pencil, the 
Casimir functions depend on the parameter A. If the rank of the Poisson pencil 
is maximal, the Casimir function is essentially unique (two Casimir functions 
are functionally dependent). The main content of the GZ theorem is that there 
exists a Casimir function depending polynomially on the parameter A, and that 
the degree of this polynomial is exactly n if dimM = 2n+l. Thus we can write 
the Casimir function in the form 

C(\) = C \ n + C 1 \ n - 1 + --- + C n . (2.15) 

This result means that the Poisson pencil selects n + 1 distinguished functions 
(Co, Ci, . . . C n ). Generically these functions are independent. Their common 
level surfaces are the leaves of the support of the pencil. Indeed, on the support 
the function C(A) must be constant independently of the particular value of A. 
Thus all the coefficients (C , C\, . . . C n ) must be separately constant. Further- 
more, as a consequence of the fact that C(A) is a Casimir function, it is easily 
seen that the coefficients C k verify the Lenard recursion relations, 

{;C k } 1 = {-,C k+1 } , (2.16) 

together with the vanishing conditions 

{•,C }o = {-,C n } 1 = 0. (2.17) 

In the language of the previous lecture, the functions (Co, C\, . . . C n ) form a 
Lenard chain. A typical property of these functions is to be mutually in invo- 
lution: 

{C j ,C k } o = {C j ,C k } 1 = . (2.18) 
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This is proved by repeatedly using the recursion relation ( |2.16|) to go back and 



forth along the chain. It follows that the leaves of the support are isotropic 
submanifolds, but since they are of maximal dimension n they are actually La- 
grangian submanifolds. These short remarks should give a sufficiently detailed 
idea of the meaning of the GZ theorem. 

Exercise 2.6 Check that that the integrals of motion Ho and H\ of the reduced 
flow X% on the invariant submanifold M3 considered in Section |1.3| are the 
coefficients of the Casimir function C(A) = XH + Hi of the Poisson pencil 
defined on M 3 . □ 



Exercise 2.7 Prove the claim (|2.18|) about the involutivity of the coefficients 



of a Casimir polynomial. □ 

From our standpoint, the above results are worthwhile of interest for two dif- 
ferent reasons: First of all they show how the Lenard recursion relations, char- 
acteristic of the theory of "soliton equations", arise in a theoretically sound 
way in the framework of bi-Hamiltonian manifold. Secondly, they highlight the 
existence of a distinguished set of Hamiltonian (Co, C\, . . . C n ) on the mani- 
fold M. Let us now choose one of the brackets of the pencil, say the bracket 
{•,-} . The function Cq is a Casimir function for this bracket, and therefore 
its level surfaces are the synplectic leaves of the bracket {-, -} . Let us call 
co>o the symplectic 2-form defined on these submanifolds. As a consequence of 
the involution relation ( [2.18|) , the restrictions of the n functions (C\, . . . C n ) 



to the symplectic leaf are in involution with respect to uj®. According to the 
Arnol'd-Liouville theorem, they define a family (or "hierarchy") of n completely 
integrable Hamiltonian vector fields on the symplectic leaf. 

Definition 2.8 The family of completely integrable Hamiltonian systems de- 
fined by the functions (Ci, . . . C n ) on each symplectic leaf of the Poisson bracket 
{■, -} will be called the GZ hierarchy associated with the Poisson pencil {-, -} A 
defined on the bi-Hamiltonian manifold M. 

We shall be particularly interested in the study of this hierarchy for two reasons. 
First we want to show that the previous simple concepts allow to reconstruct a 
great deal of the KdV hierarchy, up to the linearization process on the infinite- 
dimensional Sato Grassmannian. In other words, we want to show that the 
theory of Poisson pairs is a natural gate to the theory of infinite-dimensional 
Hamiltonian systems described by partial differential equations of evolutionary 
type. Secondly, in a finite-dimensional setting, we want to show that the GZ 
vector fields are often more than integrable in the Liouville sense. Indeed, under 
some mild additional assumptions on the Poisson pencil, they are separable, and 
the separation coordinates are dictated by the geometry of the bi-Hamiltonian 
manifold. This result strenghtens the connection between Poisson pairs and 
integrability. 
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3 A first class of examples and the reduction 
technique 



The aim of this lecture is to present a first class of nontrivial examples of GZ 
hierarchies. The examples are constructed to reproduce the reduced KdV flows 
discussed in the first lecture. The relation, however, will not be immediately 
manifest, and the reader has to wait until the fifth lecture for a full understand- 
ing of the motivations for some particular choice herewith made. 

This lecture is split into three parts. In the first one we introduce a simple 
class of bi-Hamiltonian manifolds called Lie-Poisson manifolds. They are duals 
of Lie algebras endowed with a special Poisson pencil of Lie-theoretical origin. 
The Hamiltonian vector fields defined on these manifolds admit a Lax represen- 
tation with a Lax matrix depending linearly on the parameter A. In the second 
part we show how to combine several copies of these manifolds, in such a way 
to obtain Hamiltonian vector fields admitting a Lax representation depending 
polynomially on the parameter A. Finally, in the third part, we introduce the 
geometrical technique of reduction of Marsden and Ratiu. It will allow us to 
specialize the form of the Lax matrix. The contact with the KdV theory, to 
be done in the fifth lecture, will then consist in showing that the reduced KdV 
flows admit exactly the Lax representation of the Hamiltonian vector fields con- 
sidered in this lecture. This will ascertain the bi-Hamiltonian character of the 
reduced KdV flows. The lecture ends with an example worked out in detail. 

3.1 Lie-Poisson manifolds 

In this section M = g* is the dual of a Lie algebra g. We denote by S a point 
dF 

in M, and by — — the differential of a function F : M — >• !L This differential is 

dS 

the unique element of the algebra g such that 



along any curve passing through the point S. The Poisson pencil on M is 
defined by 



where A is any fixed element in g*. In all the examples related to the KdV 

dF 

theory, g = sl(2), S and — — are traceless 2x2 matrices, and 




(3.1) 




(3.2) 







\ 



A = 



(3.3) 



V 



i o 



/ 
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The Hamiltonian vector field Xp has the form 

S 



dF 

s + " A -os 



It is already in Lax form, with a Lax matrix given by 

L(X) = XA + S . 



(3.4) 



(3.5) 



Exercise 3.1 Compute the Poisson tensor and the Hamiltonian vector fields 
associated with the pencil ( |3.2|) . □ 



3.2 Polynomial extensions 

We consider two copies of the algebra g. Accordingly, we denote by (So, Si) a 
f dF dF\ 

point in M and by I tttti 7777- the differential of the function F : M — > R. By 



v dS ' dSu 

definition, along any curve 1 1— »■ (So(t), Si(t)) we have 



dt 



, Sq\ + /tttt, S% 



The two copies of the algebra are intertwined by the Poisson brackets. As a 
Poisson pair on M we choose the following brackets 



{F,G} =(A, 
{F,G}i = (A. 



OF dG 



dS dSi 
OF dG 



+ 



OF OG 



dSi ' dSi 



dSi ' dS 

, dF <)G 

Sn 



Si, 



dF 8G 



dS n ' dS n 



dS n ' dS n 



(3.7) 



The motivations can be found for instance in |20[ (see also [0). Later on we 
shall see how to extend this definition to the case of an arbitrary number of 
copies. 

Exercise 3.2 Check that equations ( |3.7| ) indeed define a Poisson pair. □ 

Let us now study the Hamiltonian vector fields. Those defined by the brackets 
{-, -} have the form 



So 
Si 



dF 
dS~i 
dF 
dK 



+ 



Si: 



dF 
dS~ n 



(3.1 



Those defined by the second bracket {•, -}i are 

dF 



So 
Si 



So, 



dSn 



A 



dF' 



(3.9) 
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It turns out that the Hamiltonian vector fields associated with the Poisson pencil 
are given by 



Sq — — 
Si = - 



Sq + A S\ 

A A 



dF 
dS~o 









■as. 



OS, 



+ 



.4. 



OF 



(3.10) 



This computation allows to display an interesting property of these vector fields. 
If we multiply the second equation by A and add the result to the first equation 
we find 



(A 2 A + A Sx + S 



dF 
dS~o 



A 2 A + A Si + S 



(3.11) 



This is a Lax representation with Lax matrix L(X) = A 2 A-\- A Si+ Sq. It depends 
polynomially on the parameter of the pencil. We have thus ascertained that all 
the Hamiltonian vector fields relative to the Poisson pencil ( 3.1 0|) admit a Lax 
representation. The converse, however, is not necessarily true. Indeed, it must 
be noticed that the single Lax equation ( |3.11| ) is not sufficient to completely 
reconstruct the Poisson pencil ( |3.10| ). Additional constraints on the matrix 
L(A) are required to make the problem well-posed. The kind of constraints to 
be set are suggested by the geometric theory of reduction which we shall now 
outline. 



3.3 Geometric reduction 

We herewith outline a specific variant || of the reduction technique of Marsden 
and Ratiu for Poisson manifolds. This variant is particularly suitable for 
bi-Hamiltonian manifolds. 

Among the geometric objects defined by a Poisson pair (Po, Pi) on a manifold 
M we consider: 

i) a symplectic leaf S of one of the two Poisson bivectors, say Po- 

ii) the annihilator (TS)° of the tangent bundle of S, spanned by the 1-forms 

vanishing on the tangent spaces to S. 

iii) the image D = Pi(TS)° of this annihilator according to the second Poisson 

bivector P x . It is spanned by the Hamiltonian vector fields associated with 
the Casimir functions of Po by P\. 

iv) the intersection E = D n TS of the distribution D with the tangent bundle 

of the selected symplectic leaf S. 

It can be show that E is an integrable distribution as a consequence of the 
compatibility of the Poisson brackets. Therefor we can consider the space of 
leaves of this distribution, N = S/E. We assume iV to be a smooth manifold. 
By the Marsden-Ratiu theorem, iV is a reduced bi-Hamiltonian manifold. 
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The reduced brackets on N can be computed by using the process of "pro- 
longation of functions" from iV to M. Given any function / : iV — > K, we 
consider it as a function on S, invariant along the leaves of E. Then we choose 
any function F : M — > R which annihilates D and coincides with / on S. This 
function is said to be a prolongation of /. It is not unique, but this fact is 
not disturbing. It can be show that, if F and G are prolongations of / and 
g, their bracket {F, G}\ is an invariant function along E. Therefore it defines 
a function on iV which is by definition the reduced bracket {/, g}\. The final 
result, of course, is independent of the particular choices of the prolongations 
F and G. 



3.4 An explicit example 

According to the spirit of these lectures, rather than discussing the proof of the 



reduction theorem stated in Section |3.3| , we prefer to illustrate it on a concrete 
example. Let us thus perform the reduction of the Poisson pencil defined on 
two copies of q = sl(2). The matrices So and Si are traceless matrices whose 
entries we denote as follows: 



So 



( \ 

Po r 



\ 



qo -Po 



Si 



( \ 

Pi n 



\ 



Qi ~Pi 



(3.12) 



The space M has dimension six, and the entries of Sq and Si are global coor- 
dinates on it. In these coordinates the differential of a function F : M — > R is 
represented by the pair of matrices 



dF 
dS~o 



( 1 dF 
2dp~o 
OF 
dr 



OF 

dqo 
1 dF 

'2dp~o 



\ 



dF 
dS~i 



( 1 dF 
2~dp~i 

dF 

\ dri 



dF 
dqi 
1 dF 
2dp~i 



\ 



(3.13) 



/dF 

Exercise 3.3 Define the pairing \ -^,S) on q as the trace of the product 



dF 



\dS 



of the matrices and S. Show that the matrices ( |3.13| ) verify the defining 



equation (|3.6|). 



□ 
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The Hamiltonian vector fields ( p.8|) and (|3.9| ) are consequently given by 

dF dF dF 
dr dq (% 
dF dF dF 

Qo = 2 ^7T + a~ 

dpo dr dpi 

dF dF 

r = 2p 1 - n — 

9g <9p (3.14) 

dF 

P 1 = ~7T~ 
oq 

dF 

dpo 
r'i = 



and by 



dF dF 
dr dq 

dF dF 

q = 2p - q — 

dr dp 

dF dF 
r = -2po— + r — 

dq dp (3.15) 



Pi 
4l 



dF 



dqi 
dF 

dpi 

r'i = 



respectively. 

Step 1: The reduced space N. 

First we notice that the Hamiltonian vector fields ( 3.14j) verify the constraints 

r\ = 0, {r + pi + r iqi )' = . (3.16) 

It follows that the submanifold S C M defined by the equations 

n = 1 , r + p\ + nqi = (3.17) 

is a symplectic leaf of the first Poisson bivector Po- Furthermore, it follows 
that the annihilator (TS)° is spanned by the exact 1-forms dri and d(ro + p\ + 
riqi). By computing the images of these forms according to the second Poisson 
bivector ( |3.15| ), we find the distribution D. It is spanned by the single vector 
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field 



Po 


= -r 


g'o 


= 2 Po 


r'o 


= 


Pi 


= -1 


i\ 


= 2jn 




= 



(3.18) 



which verifies the five constraints 



(Po - r pi)" = 0, (g + 2p Pi - ?"oP?)' = °> 
(<?i+P?)' = 0, r = 0, ri = 0. 



(3.19) 



They show that D c TS 1 , and therefore E = D. Moreover they yield that the 
leaves of E on S are the level curves of the functions 

u x = q x + p\ 

u 2 = Po + Pigi + Pi (3.20) 
w 3 = g + 2p Pi + 9iPi + Pi 

We conclude that: 

• iV = M 3 ; 

• (mi, m 2 , M3) are global coordinates on N; 

• the canonical projection n : S — > is defined by equations (|3.20|) . 

Step 2: The reduced brackets. 

Consider any function / : N — > R. The function 

^ := fill + P?,Po +Pi<?i + P?,<?o + 2p Pi + <7iPi +Pi) (3.21) 

is clearly a prolongation of / to M, since it coincides with / on S, and is 
invariant along D. We can thus use F to compute the first component of the 
reduced Hamiltonian vector field on N according to the following algorithm: 



• (F20I) . „ . (Bill <9F <9F 
Ui^^qi + 2p 1 p 1 ^"- pi — 

9P0 9qo (3.22) 



The other components are evaluated in the same way. The final result is that 
the Hamiltonian vector fields associated with the reduced Poisson pencil on iV 
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are defined by 

ui = {ui + A)- h 2m 2 ^— 

a-u 2 ait 3 

df df 
u 2 = + A /- + (tt 3 - 2A«i)^- (3.23) 
aui 0*1/3 

w 3 = -2w 2 - h (2A«i - M3j^— 

aiti ait 2 

With this reduction process we passed from a six-dimensional manifold M to 
a three dimensional manifold N. Later on, we shall see that this manifold 
coincides with the invariant submanifold M 3 of KdV, defined by the constraint 

Uxxx — 6uu x = . (3.24) 

Step 3: the GZ hierarchy. 

To compute the Casimir function of the pencil (|3.23j) we notice that these vector 



fields obey the constraint 

(2A«i — Uz)u\ + 2w 2 ii 2 — {u\ + A)«3 = . (3.25) 

Therefore, integrating this equation, we obtain that 

C(A) = \{u\ - u 3 ) + {u\ - Ul ii 3 ) = AQ, + d (3.26) 

is the Casimir sought for. It fulfills the scheme of the GZ theorem, and it defines 
a "short" Lenard chain 

P dC = PidQ, = Podd = X 1 P 1 dC 1 = . (3.27) 

Therefore the GZ "hierarchy" consists of the single vector field 

lii = 2w 2 

Xt : u 2 = u 3 + 2ul (3.28) 
it 3 = 4iiiii 2 

As a last remark, we notice that this vector field coincides with the restriction 

du 

of the first equation — — = u x of the KdV hierarchy on the invariant sub- 

manifold ( |3.24| ). Indeed, by the procedure explained in Section [I], the reduced 
equation written in the "Cauchy data coordinates" (it, u x , u xx ) is given by 

<9tt 
du T 

= u xx (3.29) 

du xx 

— — = buu x 
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We can now pass from (|3.29|) to (|3.28| ) by the change of variables 



Mi = -u, u 2 = ~u x , u 3 = -u xx - -u 2 . (3.30) 

This remark shows that the simplest reduced KdV flow is bi-Hamiltonian. In 
the fifth lecture we shall see that this property is general, and we shall explain 
the origin of the seemingly "ad hoc" change of variables (|3.30| ). 



3.5 A more general example 



To deal with higher-order reduced KdV flows, we have to extend the class of 
bi-Hamiltonian manifolds to be considered. We outline the case of three copies 
of the algebra q. The formulas are similar to the ones of equation (|3.7|) , albeit 
a little more involved. The brackets {F, G}q and {F, G}i are now given by 



{F,G}o = (A 



" dF 


dG' 


[dS ' 


dS 2 _ 


'dF 


dG' 


_dS 


dS K 


'dF 


dG' 


[dSo 


dS \ 



+ 



dF dG 



dSi ' dSi 



+ 



dF dG 



dS 2 ' dS 



and 



{F,Gh 



A. 



+ (s. 



>2- 



Sn 



+ 

;> 


'dF 


'dF 


dG' 


_dS 1 


y dS 2 _ 


\dF 


dG' 


_dS 1 


' dSi_ 


'dF 


dG' 


[dS 


'dS \ 



(3.31) 



+ 



dF dG 



dS 2 ' as x 



(3.32) 



The comparison of the two examples allows to infer by induction the general 
rule for the Poisson pair, holding in the case of an arbitrary (finite) number of 
copies of g. The pencil ( |3.31| )- (l3.32j ) can be reduced according to the procedure 
shown before. If g = sl(2) and A is still given by ( |3.3j ), the final result of 
the process is the following: We start from a nine-dimensional manifold M 
and, after reduction, we arrive at a five-dimensional manifold N. It fulfills the 
assumption of the GZ theorem. The GZ hierarchy consists of two vector fields, 
which are the reduced KdV flows given by (|1.31|) and ( |1.32j ). 



Exercise 3.4 Perform the reduction of the pencil ( |3.31| )-( |3.32| )) for q 



sl(2). 
□ 
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4 The KdV theory revisited 



In this lecture we consider again the KdV theory, but from a new point of view. 
Our purpose is twofold. The first aim is to show that the KdV hierarchy is an- 
other example of GZ hierarchy. The second aim is to explain in which sense the 
KdV hierarchy can be linearized. The algebraic linearization procedure dealt 
with in this lecture was suggested for the first time by Sato |27| (see also the 
developments contained in |7|, ^8|, p9[), who exploited the so-called Lax repre- 



sentation of the KdV hierarchy in the algebra of pseudo-differential operators. 
Here we shall give a different description, strictly related to the Hamiltonian 
representation of the KdV hierarchy as a kind of infinite-dimensional GZ hi- 
erarchy. However, the presentation does not go beyond the limits of a simple 
sketch of the theory. We refer to [110] for full details. 



4.1 Poisson pairs on a loop algebra 



In this section we consider the infinite-dimensional Lie algebra M of C°°-maps 
from the circle S 1 into q = sl(2). A generic point of this manifold is presently 
a 2 x 2 traceless matrix 



S 



^ p(x) r(x) ^ 



q(x) —p(x) 



(4.1) 



/ 



whose entries are periodic functions of the coordinate x running over the circle. 
The three functions (p, q, r) play the role of "coordinates" on our manifold. The 
scalar-valued functions F : M — > R to be considered are local functionals 



f(p,q,r;p x ,q x ,r x ;...)dx 



(4.2) 
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As before, their differentials are given by the matrices 

\ 



5F 

Is 



( l Jl 

2 5p 

Sf 
\ 5r 



Sf 
6q 

2Sp / 



(4.3) 



whose entries are the variational derivatives of the Lagrangian density / with 
respect to the functions (p,q,r). The Poisson pencil is similar to the first one 
considered in the previous lecture (see equation (|3.2j)). It is defined by 



{F,G}> 



S + XA, 



6£ 6G 
IS' 5S 



5F_ 6G 

1 Is' Is 



(4.4) 



It differs from the previous example by the addition of the nontrivial cocycle 

da 



u (a, b) 



1 dx 



■bdx 



(4.5) 
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This term is essential to generate partial differential equations. It is responsible 
for the appearance of the partial derivative in the expansion of the Hamiltonian 
vector fields 



S= ^ 



5F_ 

Is 



S 



(4.6) 



Exercise 4.1 Recall that a two-cocycle on g is a bilinear skewsymmetric map 
w:gxg^R which verifies the cyclic condition 

u(a, [b, c]) + u(b, [c, a]) + lu(c, [a, b}) = . 

Using this identity and the periodic boundary conditions check that equation 
defines a Poisson bivector. □ 



4.2 Poisson reduction 

We apply the same reduction technique used in the previous lecture, avoiding 
to give all the details of the computations. They can be either worked out by 
exercise or found in ||, 19 



The first Poisson bivector Pq is defined by 



S 



A ' SS 



(4.7) 



where A is still defined by equation (|3.3| ). These Hamiltonian vector fields obey 
the only constraint r = 0. Therefore the submanifold S formed by the matrices 



S 



\ 



-P 



(4.8) 



/ 



is a symplectic leaf of Pq. The annihilator (TS)° is spanned by the differentials 
of the functionals F : M — > R depending only on the coordinate function r. 
Consequently, the distribution D is spanned by the vector fields 



P 



Sf 
8r 



5r) x Sr 



(4.9) 







The distribution D is thus tangent to S and E coincides with D. The vector 
field (14.91) verifies the constraint 



q + 2pp + p x = (q + p 2 + p x y = . 



(4.10) 
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It follows that the leaves of the distribution E are the level sets of the function 



u = q + p + p a 



(4.11) 



Therefore the quotient space N is the space of scalar functions u : S 1 — > K, and 
(|4.11|) is the canonical projection 7r : S — > S/E. We see that the manifold N is 
(isomorphic to) the phase space of the KdV equation. 

We use the projection ( |4.11| ) to compute the reduced Poisson bivectors. The 
scheme of the computation is always the same. First we prolong any functional 
T = f sl f(u, u x , ■ ■ ■ )dx on N into the functional 



F(p, q, r) 



f(q + p 2 + p x , q x + 2pPx + Pi 



on S. Then we compute its differential at the points of S, 



5_F_ 

5S 



( I /6F\ SF_ 
~2 1^7 J +P fa 



5F_ 
Su 



dx 



\ 



V 



o 



1 (8F\ SF 

2 \J^) X ~ P 1^ ) 



(4.12) 



(4.13) 



Finally, we evaluate the reduced Hamiltonian vector fields on iV according to 
the usual scheme: 



.(4-11). . „ . 

u = q + p x + 2pp 



(4.6) 



Sf 
5u 



+ 9 + A)— -2p— 
op or 



+ 2p 

(4.13) 



+ 



(q + A) 



8f\ +2p S_f 
5u J „ 8u 



+ 



\( 5 1 
2 \5p 



lf 5 l 
2\Su 



+ j- + (q + X)j- 
or oq 



o~ f 

+ {qx + p X x + 2pp x )— 

OU 



(4.14) 



We obtain the Poisson pencil of the KdV equation. This pencil is therefore the 
reduction of the "canonical" pencil (|4.4j) over a loop algebra. 



4.3 The GZ hierarchy 

The simplest way for computing the Casimir function of the above pencil is to 
use the Miura map. Since this map relates the pencil to the simple bivector 
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of the mKdV equation, it is sufficient to compute the Casimir of the latter 
bivector, and to transform it back to the phase space of the KdV equation. 
We notice that the Casimir function of the mKdV hierarchy (|1.39[) is given 

by 

H(/i) = 2z / hdx , (4.15) 
Js 1 

where the constant z has been inserted for future convenience. 

To obtain the Casimir function of the KdV equation, we must "invert" the 
Miura map by expressing h as a function of u. To do that we exploit the 
dependence of the Miura map on the parameter A = z 2 of the pencil. We 
know that in the finite-dimensional case the Casimir function can be found 
as a polynomial in A. In the infinite-dimensional case, we expect the Casimir 
function to be represented by a series. It is then natural to look at h in the 
form of a Laurent series in z, 

h(z) = z + J2 h i z ~ l > ( 4 - 16 ) 

i>i 

whose coefficients hi are scalar-valued periodic functions of x. In this way we 
change our point of view on the Miura map. Henceforth it must be looked at as 
a relation between a scalar function u and a Laurent series h(z). This change 
of perspective deeply influences all the mKdV theory. It is a possible starting 
point for the Sato picture of the KdV theory, as we shall show later. 

By inserting the expansion (|4.16| ) into the Miura map h x + h 2 = u + z 2 and 



equating the coefficients of different powers of z, we easily compute recursively 
the coefficients hi as differential polynomial of the function u. The first ones 
are 



h = 


\u 


h 2 = 


~i U x 


h = 


g i^xx U ) 


hi = 




h 5 = 


"^i^xxxx QuU X x 



(4.17) 



5u 2 x + 2u 3 ) . 

One can notice (see [l[]) that all the even coefficients h 2 i are total x-derivatives. 
This remark explains the "strange" enumeration with odd times used for the 
KdV hierarchy in the first lecture. 

To compute concretely the GZ vector fields, besides the Casimir function 



H{u,z) = 2z^2l hiz- l dx, (4.18) 
i>i J sl 
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we need its differential. To simplify the notation we set 



a:= — = l + Ya lZ - 1 . (4.19) 

ou ^ 

v>\ 

Once again, the simplest way for evaluating this series is to use the Miura map. 
We notice that j3 = 2z is the differential of the Casimir of the mKdV equation. 
From the transformation law of 1-forms, 

*ft» = /9, ( 4 - 2 °) 
we then conclude that that a solves the equation 

-a x + 2ah = 2z . (4.21) 

As before, the coefficients a\ can be computed recursively. One finds a Laurent 
series in A = z 2 , 

a = 1 - 1-uX- 1 + ~ (3m 2 - u xx )\~ 2 + ■■■ , (4.22) 

/ o 

whose first coefficients have already appeared in (|1.20|) . From a we can easily 
evaluate the Lenard partial sums ay> = (A J a) + and write the odd GZ equations 
in the form 



du 



m 2j . 

The even ones are 

du 



- = ( - \d xxx + 2{u + \)d x + u x ) ( Q W) . (4.23) 



dt 



2j 



. (4.24) 



The above equations completely and tersely define the KdV hierarchy from the 
standpoint of the method of Poisson pairs. 



4.4 The Central System 

We shall now pursue a little further the far-reaching consequences of the change 
of point of view introduced in the previous subsection. According to this new 
point of view, the mKdV hierarchy is defined in the space C of the Laurent 
series in z truncated form above. This affects the whole picture. 

Let us consider again the basic formulas of the mKdV theory. They are the 
Miura map, 

h x + h 2 =u + z 2 , (4.25) 
the formula for the currents ( |1.40| ), 

H {2]+1) = -ia£'> + a (]) h, H W) = , (4.26) 
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and the definition of the mKdV hierarchy 



(4.27) 



They were obtained in the first lecture. Presently they are complemented by 
the information that h(z) is a Laurent series of the form Q4.16I) . We shall now 
investigate the meaning of the above formulas in this new setting. 

We start form the series h(z), and we associate with it a new family of 
Laurent series W\z) defined recursively by 



h U+i) = h (j) + hh U) 



(4.28) 



starting from = 1. They form a moving frame associated with the point 
h in the space of (truncated) Laurent series. The first three elements of this 
frame are explicitly given by 



= h, =h x + h 2 



(4.29) 



We see the basic block h x + h 2 of the Miura transformation appearing. We call 
H + the linear span of the series {h^'}j>o. It is a linear subspace of C, attached 
to the point h. We can now interpret the three basic formulas of the mKdV 
theory as properties of this linear space: 



The Miura map ( [4.25 ) tells us that the linear space 7i + is invariant with 
respect to the multiplication by A, 



\{n+) c n. 



(4.30) 



The formula ( 4.26[ ) for the currents then entails that the currents for 
j e N, belong to H + : 



H {3) e H- 



(4.31) 



Furthermore, in conjunction with equation ( 4.21|) , it entails that the asymp- 
totic expansion of the currents has the form 



z j + 0(z- 1 ) . 



(4.32) 



i>i 



Finally, the mKdV equations ( |4.27| ) can be seen as the commutativity 

d 

conditions of the operators (d x + h) and (-^- + H^): 



d x + h,^- + H® 



(4.33) 



32 



Q 

Used together, conditions ( |4.31| ) and ( |4.33| ) imply that the operators (— — \-H^ 

Olj 

leave the linear space 7i + invariant: 

{JL+H®){H + )<ZH+ ■ (4.34) 

This is the abstract but simple form of the laws governing the time evolution of 
the currents These equations are the "top" of the KdV theory, and form 

the basis of the Sato theory. It is not difficult to give them a concrete form. By 
using the form of the expansion ( [4.32 ) it is easy to show that equations ( |4.34 ) 



are equivalent to the infinite system of Riccati-type equations on the currents 
//'': 

Q TT(j) k 

<m_ + H (j) H (k) = H U + k) + C*"0 + H(H^ . (4.35) 

h i=i i=i 

It will be called the Central System (CS). 

Exercise 4.2 Prove formulas fl4.31|) and ( |4.32j ). 



4.5 The linearization process 

The first reward of the previous work is the discovery of a linearization process. 
The equations ft4-35|) of the Central System are not directly linearizable, but they 
can be easily transformed into a new system of linearizable Riccati equations 
by a transformation in the space of currents. This idea is realized once again 
by a "Miura map" . The novelty, however, is that this map is now operating on 
the space of currents rather than on the phase space of the KdV equation. 

We simply give the final result. Let us consider a new family of currents 
{W^} k>o of the form 

W ik) = z k + J2 Wfz- 1 , (4.36) 
i>i 

and let us denote by W+ their linear span in C We define (see also p9| ) a new 
system of equations on the currents by imposing the "constraints" 

(JL + z *)( W+ )cW + (4.37) 

Otk 

on their linear span W+. It is easily seen that they take the explicit form 

^ — + z j W {k) = W i3+k) + J2 W i kwij ~ l) ■ ( 4 - 38 ) 
dt i i=i 

They will be called the Sato equations (on the "big cell of the Sato Grassman- 
nian"). They are a system of linearizable Riccati equations. This can be seen 
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either from the geometry of a suitable group action on the Grassmannian [JF| or 
by means of the following more elementary considerations. We write equations 
(|4.38|) in the matrix form 



<9W 



w • T A k - A k ■ w = wr fe w , 



(4.39) 



where W = (Wj fc ) is the matrix of the components of the currents W^ k \ A is 
the infinite shift matrix 



A 



1 ••• 
1 



(4.40) 



and T h is the convolution matrix of level k, 

••• 10 

1 1 ••• 

1 



(4.41) 



One can thus check that the matrix Riccati equation ( 4.39j ) is solved by the 
matrix 



W = V- IT 1 , 

where U and V satisfy the constant coefficients linear system 



d_ 
dt k 



u = T A fc u - r fe v 



d_ 



V = A fc V 



(4.42) 



(4.43) 



The closing remark is that the Sato equations are mapped into the Central 
System ( |4.35[) by the following algebraic Miura map: 



WW 



(4.44) 



The outcome of this long chain of extensions and transformations is the following 
algorithm for solving the KdV equation: 
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i) First we solve the linear system ( 4.43 ), with a suitably chosen initial condi- 
tion, which we do not discuss here; 



ii) Then we use the projective transformation (|4.42|) and the Miura map (|4.44| ) 

to recover the currents 

iii) Finally, we extract the first current H^ 1 ' = h, and we evaluate the first 

component hi of its Laurent expansion in powers of z~ l . 



The function 

u{x,t 3 ,...) = 2/ii| tl=JC (4.45) 
is then a solution of the KdV equation. 



4.6 The relation with the Sato approach 



The equations (|4.27| ) make sense for an arbitrary Laurent series h of the form 



( |4.16| ), even if it is not a solution of the Riccati equation h x + h = u + z . 
Hence they define, for every j, a system of PDEs for the coefficients hi. We 
will show[] that these systems are equivalent to the celebrated KP hierarchy of 
the Kyoto school (see the lectures by Satsuma in these volume). The usual 
definition of the KP equations can be summarized as follows. Let ^>T> be the 
ring of pseudodifferential operators on the circle. It contains as a subring the 
space V of purely differential operators. Let us denote with (•), the natural 
projection from ^>T> onto T>. Let Q be a monic operators of degree 1, 



Q = d 



-j 



(4.46) 



The KP hierarchy is the set of Lax equations for Q 

^-Q = [(Q j ) + ,Q\- 



(4.47) 



The aim of this subsection is to show that such a Lax representation just arises 



as a kind of a Euler form of the equations (|4.27| ). Before stating the next result, 
we must observe that the relations ( |4.28| ) can be solved backwards, in such a 
way to define the Faa di Bruno elements h^' for all j G Z. 



Proposition 4.3 Suppose the series h of the form ( \4-lo\ ) to evolve according 
to a conservation law, 



dh 
dt 



d x H, 



(4.48) 



'See also the papers ffl, t52 
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for an arbitrary H. Then the Fad di Bruno elements hP\ for j 6 7L, evolve 
according to 

(§r + H)h^^(i)(d k x H)h^ , (4.49) 

where 

'j\ j(j ~ 1) • • • (j - k + 1) fj\ 



kj k\ \0 

Now we consider the map <p : £ ^ $?T>, from the space of Laurent series to 
the ring of pseudodifferential operators on the circle, which acts on the Faa di 
Bruno basis according to 

0(7i(i)) = d j . (4.50) 

This map is then extended by linearity (with respect to multiplication by a 
function of x) to the whole space C 

Definition 4.4 We call Lax operator of the KP theory the image 

Q = (j>{z) (4.51) 
of the first element of the standard basis in C. 

If the qj are the components of the expansion of z on the Faa di Bruno basis, 

z = h {1) - , (4.52) 



then we can write 



Q = d-J2 (li d ' j ( 4 - 53 ) 

3>1 



according to the definition of the map 0. We note that equation ( |4.52j ) uniquely 



defines the coefficients qj as differential polynomials of the components hj of 
h(z): 

qi = h, q 2 = h 2 , q 3 = h + h\ 

q 4 = h 4 + 3hh 2 - hh lx (4.54) 



This is an invertible relation between the hj and the qj, so that equation (|4.52j ) 
may be seen as a change of coordinates in the space C. 

Proposition 4.5 The map (j) has the following three properties: 
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i) Multiplying a vector of the Fad di Bruno basis by a power z k of z yields 

<f)(z k ■ h {j) ) = d j ■ Q k . (4.55) 

ii) The evolution along a conservation law of the form 



translates into 

^ i {<P(h^))=Y J [d\H k }.Q k . (4.56) 

k 

iii) // 7r + and IT + are respectively the projection onto the positive part 7i + C C 
and V C ^V, then 

<p o 7T+ = n+ O <j) . (4.57) 

To obtain the Sato form of the equations (|4.27|) , we derive the equation 

z = - U h( ~ l) ( 4 - 58 ) 
i>i 

with respect to the time tj, getting 

i>i 3 3 i>i ■> 

Applying the map to both sides of this equation we obtain 

Z>1 fc>l fc>l 

or 

|^ + ^[g,^]Q- fc = 0. (4.61) 
fe>i 

Finally, we introduce the operator 

= = ( ^ + ^ j = Q 3 + ^ ^g- fc (4.62) 

V fc>i / fc>i 

associated with the current density and we note that 

flW = </>(tt + (^))) = (0(^')) + = (#) + . (4.63) 
Thus we can write ( |4.61| ) in the final form 



ff + [^(^)J = > ( 4 - 64 ) 



which coincides with equation ( |4.47| ). 
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5 Lax representation of the reduced KdV flows 



In this lecture we want to investigate more accurately the properties of the 
stationary KdV flows, that is, of the equations induced by the KdV hierarchy 
on the finite-dimensional invariant submanifolds of the singular points of any 
equation of the hierarchy. Examples of these reductions have already been 
discussed in the first lecture. In the third lecture we realized, in a couple of 
examples, that the reduced flows were still bi-Hamiltonian. Although not at 
all surprising, this property is somewhat mysterious, since it is not yet well 
understood how the Poisson pairs of the reductions are related to the original 
Poisson pairs of the KdV equation. Moreover, even if the subject is quite 
old and classical (see, e.g., §, it was still lacking in the literature a 
systematic and coordinate free proof that such reduced flows are bi-Hamiltonian 
(see, however, ]2|, [SO)). In this lecture we will not provide such a proof, which is 
contained in [11], but we will give a sufficiently systematic algorithm to compute 
the reduced Poisson pair. This algorithm is based on the study of the Lax 
representation of the reduced equations. 



5.1 Lax representation 

In this section we associate a Lax matrix (polynomially depending on A) with 
each element This matrix naturally arises from a change of basis in the 

linear space 7i + attached to the point h. So far we have introduced two bases: 

i) The moving frame {h^}; 

ii) The canonical basis {H^}. 

Presently we introduce a third basis by exploiting the constraint 

\(n+) c n + , (5.i) 

characteristic of the KdV theory. The new basis is formed by the multiples 
{X^H^, X>H^} of the first two currents. Formally we define 

iii) the Lax basis: (X J , X J h) . 

The use of this basis leads to a new representation of the currents where 
each current is written as a linear combination of the first two, = 1 and 
= h, with coefficients that are polynomials in A. Let us consider a few 
examples: 

= l + 0-h 
H {1) = • 1 + 1 • h 
#( 2 ) = A • 1 + • h 

(5 2) 

H {3) = -h 2 -l + (X-h 1 )-h 
= A 2 • 1 

#( 5 ) = (-Xh 2 + h x h 2 - h 4 ) ■ 1 + (A 2 - Xhi + h\-h z )-h . 
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This new representation also affects our way of writing the action of the oper- 
d 

ators (— h H^'). Let these operators act on and For the basic 

invariance condition ( [4.34[ ), we get an element in 7i + which can be represented 
on the Lax basis. As a result we can write 



1 




1 






h 




h 



(5.3) 



where L^\X) is the Lax matrix associated with the current We shall see 

below the explicit form of some of these matrices. 

It becomes now very easy to rewrite the Central System in the form of 
equations on the Lax matrices L^'(X). We simply have to notice that the 
equations QL35|) entail the "exactness condition" 



dt k 



(5.4) 



d d 

from which it follows that the operators f— — \-H^A and (— \-H^ k ^ commute: 

utj utk 



dtj ' dt k 







(5.5) 



It is now sufficient to evaluate this condition on (H^°\ H^) and to expand on 
the Lax basis to find the "zero curvature representation" of the KdV hierarchy: 



dti 



+ =0 . 



(5.6) 



Suppose now that we are on the invariant submanifold formed by the singular 
points of the j—th member of the KdV hierarchy. On this submanifold 







(5.7) 



and the zero curvature representation becomes the Lax representation 



[LM,L®] . 



(5i 



We have thus shown that all the stationary reductions of the KdV hierarchy 
admit a Lax representation. As a matter of fact, this Lax representation coin- 



cides JTI] with the Lax representation of the GZ systems on Lie-Poisson man- 
ifolds studied in Section 3. The latter are bi-Hamiltonian systems. Therefore, 
we end up stating that the stationary reductions of the KdV theory are bi-Ham- 
iltonian, and we can construct the associated Poisson pairs. We shall now see 
a couple of examples. 
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5.2 First example 



We study anew the simplest invariant submanifold of the KdV hierarchy, defined 
by the equation 

u X xx — Quu x = . (5-9) 

In this example we consider the constraint from the point of view of the Central 
System. Since the constraint is the stationarity of the time t 3 , we have to 
consider only the first three Lax matrices. As for the matrix L^ 1 ', the following 
computation, 

d 

+ # (1) )l = 0-l + l-/i 



d 



(» + hV)H^= ] H® + 2h 1 { =\\ + 2h) -1 + 0-h 

uti 



(5.10) 



shows that 



1 

A + 2/ix 



(5.11) 



Similarly, the computation 
d 



(^ + H^)l i ^ ) -h 2 .l + (X-h 1 ).h 
ot 3 

(JL + H m H d)(^) H (A) + hiHi 2) + h2H d) + h + H s 

y dt 3 

{ =\x 2 + A/i! + 2h 3 - h\) • 1 + h 2 ■ h 



(5.12) 



yields 



L (3) 



V 



(5.13) 



/ 



—h 2 X — hi 

X 2 + hX + 2h 3 - h 2 h 2 

On the submanifold M 3 defined by equation ( |5.9|) this matrix verifies the Lax 
equation 



[LM,L®] . 



(5.14) 



This equation completely defines the time evolution of the first three compo- 
nents (hi,h 2 ,h 3 ) of the current = h. These components play the role of 
coordinates on M 3 . We get 



dh 
dh 
dti 
dh 
dh 



-2h 2 

-2h 3 - h 2 
-Ahih 2 



(5.15) 
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By the change of coordinates 



hi 



-u 



ho = — u a 
4 



coming from the inversion ( [4.15 ) of the Miura map, these equations take the 
form 



du 



du x 



du : , 



6uu n 



(5.16) 



already encountered in Lecture [IJ. This shows explicitly the connection between 
the two points of view. 

To find the connection between these equations and the GZ equations dealt 
with in the first example of Lecture ^ we compare the Lax matrix 



L' 3 >(A) = A ! 



1 



+ A 



1 

h x 



-h 2 -hi 

,2 



^ 2h 3 — h{ h 2 



with the Lax matrix 



S{\) = A 2 



(° 


0^ 


+ A 


/ 

Pi 


1 ^ 


+ 


( 

Po 


-(<li+P 2 i) 


\ 


V 


°y 






-Pi J 




y % 


-Po 


J 



associated with the points of the symplectic leaf defined by ( |3.17| ). We easily 
identify L^> with the restriction of S(X) to pi = upon setting 



Po 



-h 2 , qi = h, q = 2h 3 - h\ . 



(5.17) 



By comparing these equations with the projection Q3.20| ), which allows to pass 
from the symplectic leaf S to the quotient space N = S/E, we obtain the change 
of coordinates 



ui = h h 



u 2 



-h 



2- 



u 3 = 2h 3 - hi, 



(5.18) 



connecting the reduction ( |5.15| ) of the Central System to the GZ system ( |3.28j ) 
dealt with in the third Lecture. The latter was, by construction, a bi-Hamil- 
tonian system. We argue that also the reduction of the Central System herewith 
considered is a bi-Hamiltonian vector field, and that its Poisson pair is obtained 
by geometric reduction. Basic for this identification is the property of the Lax 
matrix of being a section of the fiber bundle 7r : S — > S/E appearing in the 
geometric reduction. It is this property which allows to set an invertible relation 
among the coordinates {ui,u 2 ,u 3 ), coming from the geometric reduction, and 
the coordinates [hi, h 2 , h 3 ) coming from the reduction of the Central System. 
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5.3 The generic stationary submanifold 

It is now not hard to give the general form of the matrices for an arbitrary 
odd integer 2j + 1. First we observe that 

d 



(-: 



Then we notice that 
d 



(5.19) 



ot ■ / — ' 



2 v * 



(5.20) 



i=i 



Therefore 



( 



a 



+ ELiWr 1) ) + # 1 a^ +1 ) + EL^« ' _1) y 



(5.21) 



By using the definition ( |4.21 ) of the Lenard series a(z) of which the polynomials 
av) are the partial sums, it is easy to prove that is a traceless matrix. 

We leave to the reader to specialize the matrix and to write explicitly 
the Lax equations 



dtx 



[LW,L®] 



dU 



(5.22) 



They should be compared with the reduced KdV equations ( 1.31|) and ( |1.32j ) 
on the invariant submanifold defined by the constraint 



— 10uu xxx — 20u x u xx + 30m u a 







(5.23) 



They should also be compared with the GZ equations obtained via the geometric 
reduction process applied to the Lie-Poisson pairs defined on three copies of 
st(2) by equations (|3.31|) and (|3.32|) . We have not displayed explicitly these 
equations yet. We will give their form in the next lecture. 



5.4 What more? 

There is nothing "sacred" with the KdV theory. As we know, it is related with 
the constraint 



z 2 (H+) C (H 



+) , 



(5.24) 



which defines an invariant submanifold of the Central System. Many other 
constraints can be considered. For instance, the constraint 



z 3 (H + ) c (H+) 



(5.25) 
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leads to the so-called Boussinesq theory, and is studied in fl^] . What is re- 
markable is that the change of constraint does not affect the algorithm for the 
study of the reduced equations. All the previous reasonings are valid without 
almost no change. The only difference resides in the fact that the computations 
become more involved. This remark allows to better appreciate the meaning of 
the process leading from the KdV equation to the Central System. We have not 
only given a new formulation to known equations. We have actually found a 
much bigger hierarchy, possessing remarkable properties, which coincides with 
the KdV hierarchy on a (small) proper invariant subset. The integrability prop- 
erties belong to the bigger hierarchy, and hold outside the KdV submanifold. 
Many other interesting equations can be found by other processes of reduction. 
There is some evidence that a very large class of evolution equations possessing 
some integrability properties can be eventually recovered as a suitable reduc- 
tion of the Central System, or of strictly related systems. However, we shall not 
pursue this point of view further, since it would lead us too far away from our 
next topic, the separability of the reduced KdV flows. 
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6 Darboux— Nijenhuis coordinates and Separa- 
bility 

In this lecture we shall consider the reduced KdV flows from a different point of 
view. Our aim is to probe the study of the geometry of the Poisson pair which, 
as realized in the third and fifth lectures, is associated with these flows. The 
final goal is to show the existence of a suitable set of coordinates defined by and 
adapted to the Poisson pair. They are called Darboux-Nijenhuis coordinates. 
We shall prove that they are separation coordinates for the Hamilton- Jacobi 
equations associated with the reduced flows. 

To keep the presentation within a reasonable size, we shall mainly deal 
with a particular example, and we shall not discuss thoroughly the theoretical 
background, referring to ]TTJ for more details. We shall use the example to 



display the characteristic features of the geometry of the reduced manifolds. The 
reader is asked to believe that all that will be shown is general inside the class 
of the reduced stationary KdV manifolds, whose Poisson pencils are of maximal 
rank. A certain care must be used in trying to extend these conclusions to other 
examples like the Boussinesq stationary reductions, whose Poisson pencils are 
not of maximal rank. They will not be covered in these lecture notes. The 
example worked out is the reduction of the first and the third KdV equations 
on the invariant submanifold defined by the equation 



10uu xxx — 20u x u xx + 30u 2 u x = , (6.1) 



a problem addressed at the end of Section 5.3. 



6.1 The Poisson pair 

As we mentioned several times, the invariant submanifold M 5 defined by equa- 
tion (|6.1| ) has dimension five. From the standpoint of the Central System, it is 
characterized by the two equations 

5 

z 2 (H+) C (H+) , H^h = X 3 + J2 hiH^~ l) + H\ . (6.2) 

i=i 

We recall that the first constraint means that, inside the big cell of the Sato 
Grassmannian, we are working on the special submanifold corresponding to 
the KdV theory. The second constraint means that, inside the phase space of 
the KdV theory, we are working on the set (|6.1|) of singular points of the fifth 
flow. The two constraints play the following roles. The first constraint sets 
up a relation among the currents H^>: All the currents are expressed as linear 
combinations (with polynomial coefficients) of the first two currents H^ = 1 
and H^ = h. So this constraint drastically reduces the number of the unknowns 
Hf to the coefficients hi of h. The second constraint then further cuts the degrees 
of freedom to a finite number, by setting relations among the coefficients hi. It 
can be shown that only the first five coefficients (hi, h 2 , h 3 , hi, h 5 ) survive as free 
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parameters. All the other coefficients can be expressed as polynomial functions 
of the previous ones. By a process of elimination of the exceeding coordinate, 
one proves that the restriction of the first and third flows of the KdV hierarchy 
are represented by the following differential equations: 



and 



3 



dh 



at 



3 



at. 



-2ho 



dhx 
dh 

— = -2h 3 - hi 
at % 

dh 3 



dh 



-2h x h 2 - 2/i 4 (6.3) 



™± = -2h 5 - h 2 2 - 2h 1 k s 
at % 

Oils, n 

— i = -Ah 3 h 2 + 2h 2 h 2 - Ah x h A 
ati 



-2h 4 + 2h 1 h 2 



2 - -2h 5 + h 2 2 + 



-2/ii/i 4 + tti x l h 2 - 2h 3 h 2 (6.4) 



dh 

= -2h 3 2 - 2h 2 h A + 2h x h 2 2 + + h x 2 h z 

dt 3 
dh 5 



2h 2 h A - Ah 3 h A + 2h^h 2 



They can also be seen as the Lax equations ( |5.22j ). However, for our purposes, 
it is more important to recognize that the above equations are the GZ equations 
of the Poisson pencil defined on M 5 . This pencil can be computed according 
to the reduction procedure explained in the third lecture. The final outcome is 
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that the reduced Poisson bivector is given by 
dH „ ^dH , dH 

<9if „ , ,<9F 3# 9x <9# 



dH dH dH 

h 3 = 2(2/n - A)— + (2/13 + 2/i? - 4A/ il )— + 2(/i 4 + 

cm 2 CW4 <%5 

^4 = 2(A ~hi)-^ + 2h 2 — - (2h 3 + 2h\ - 4A/H)— (6.5) 

+ (2h 5 - 6/li/lg + ^2 + 2/ *i + 4A/i 3 + 2A/i?) — 

a/15 

h = ~2h 2 -^- + (2h 3 + h\- ^h x )— - 2(/i 4 + h M-Q^ 

Si 1 / 

- (2h 5 - Qhths + h 2 2 + 2h\ + 4A/i 3 + 2A/i 2 ) — . 

C//14 

The Casimir function of this pencil is a quadratic polynomial, 

C(A) = Q,A 2 + C X A + C 2 , (6.6) 

and the coefficients are 
C = - 2/11/13 + h s 

Ci = /i 2 /i 4 - M5 + |/ii 2 ^3 - \h x h 2 - \h 3 - \h x 
C 2 = \h 3 h 2 - h 3 h 5 + i/i! 5 + h x h 3 - h x h 2 h A - ^h 3 h 3 + /ii 2 /i 6 + \h 



2 



(6.7) 



The Lenard chain is 



P dC = 

dh 



P dC 2 = PidCi = ^ 

dt 3 



^ (6. 



PidC 2 = 



where h is the vector (hi, h 2 , h 3 , h±, hs). It shows that the reduced flows are 
bi-Hamiltonian. Finally, if one uses the coordinate change ( [4.171 ) from the co- 
ordinates (hi, h 2 , h 3 , hi, h 5 ) to the coordinates (u,u x ,u xx ,u xxx ,u xxxx ), one can 
put the equations (|6.3|) and (|6.4j ) in the form (|1.31| ) and ( |1.32|) considered in 
the first lecture. 
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6.2 Passing to a symplectic leaf 



We aim to solve equations ( |6.3| ) and ( |6.4| ) by the Hamilton-Jacobi method. This 
requires to set the study of such equations on a symplectic manifold. This can 
be easily accomplished by noticing that these vector fields are already tangent 
to the submanifold S4 defined by the equation 

Co = E , (6.9) 

for a constant E. We know that this submanifold is symplectic since Cq is the 
Casimir of P . The dimension of S4 is four, and the variables (hi, h 2 , h 3 , h^) 
play the role of coordinates on it. 

For our purposes it is crucial to remark an additional property of S4: It is a 
bi-Hamiltonian manifold. This means that also the second bivector Pi induces, 
by a process of reduction, a Poisson structure on S4 compatible with the natural 
restriction of Pq. This is not a general situation. It holds as a consequence of a 
peculiarity of the Poisson pencil (|6.5| ). The property we are mentioning concerns 
the vector field 

Z=A ( 6 .io) 

One can easily check that: 

i) Z is transversal to the symplectic leaf £4. 

ii) The functions which are invariant along Z form a Poisson subalgebra with 

respect to the pencil. 

In simpler terms, the Poisson bracket of functions which are independent of h§ 
is independent on h 5 as well. Since they coincide with the functions on S 4 (by 
the transversality condition), this property allows us to define a pair of Poisson 
brackets also on S^. The first bracket is associated with the symplectic 2-form 
u>o on S4. It can be easily checked that 

ujq = h\dh\ A dh,2 + —{dh-i A dh^ + dh§ A dhi) . (6-11) 

The second Poisson bracket can be represented in the form 

{f,g} 1 = u; (NX f ,X g ) , (6.12) 

where Xf and X g are the Hamiltonian vector fields associated with the functions 
/ and g by the symplectic 2-form u , and N is a (1, l)-tensor field on S4, called 
the Nijenhuis tensor associated with the pencil (see, e.g., [0). In our example 
one obtains 

N =( - "^ - + ^ - ~ 2h ^w) 8 

+ (*. _ hi) ± dh 2 + (A + 2 hl JL + h2 A) 9 dh3 (6 . 13) 
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Thus we arrive at the following picture of the GZ hierarchy considered in this 
lecture. It is formed by a pair of vector fields, X\ and X 3 , defined by ( |6.3| ) 
and Q6. 4| )- They are tangent to the symplectic leaf (S±,uo) defined by equations 
(|6.9|) and ( p. 11 ). This symplectic manifold is still bi-Hamiltonian, and therefore 
there exists a Nijenhuis tensor field N, defined by equation (|6.12 ). The vector 
fields Xi and X 3 span a Lagrangian subspace which is invariant with respect to 
N. One finds that they obey the following "modified Lenard recursion relations" 



NX, = X 3 + 
NX 3 = + 

From them we can extract the matrix 

/ 



(Tr N)Xt 
-det N)Xi . 



(6.14) 



V 



Tr N 1 
-det N 



/ 



(6.15) 



which represents the action of N on the abovementioned Lagrangian subspace. 
It will play a fundamental role in the upcoming discussion of the separability of 
the vector fields. 

Exercise 6.1 Compute the expression of the reduced pencil on S4 and check 
the form of the Nijenhuis tensor, as well as the modified Lenard recursion rela- 
tions flBmp . 



6.3 Darboux— Nijenhuis coordinates 

We are now in a position to introduce the basic tool of the theory of sepa- 
rability in the bi-Hamiltonian framework: The concept of Darboux-Nijenhuis 
coordinates on a symplectic bi-Hamiltonian manifold, like S4. 

Given a symplectic 2-form uq and a Nijenhuis tensor N coming from a Pois- 
son pencil defined on a 2n-dimensional manifold A4, under the assumption that 



the eigenvalues of N are real and functionally independent, one proves |L8| the 
existence of a system of coordinates (Ai, . . . , A n ; /ii, ... , \i n 
for u , 



which are canonical 



u 



i=l 



and which allows to put N* (the adjoint of N) in diagonal form: 



N*dXi = \4\ 



N*dfii = Xidfii 



(6.16) 



(6.17) 



The coordinates Aj are the eigenvalues of N*, and therefore can be computed 
as the zeroes of the minimal polynomial of N: 



\ n + dA"' 1 + • • • + Cn = 



(6.18) 
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The coordinates fij can be computed as the values that a conjugate polynomial 

/WiA"" 1 + ••• + /„ (6.19) 
assumes on the eigenvalues Xj, that is, 

^■ = /iAr 1 + - •• + /»> J = l,---,n. (6.20) 

The determination of this polynomial, which is not uniquely defined by the 
geometric structures present in the theory, requires a certain care. Although 
there is presently a sufficiently developed theory on the Darboux-Nijenhuis 
coordinates and on their computation, for the sake of brevity we shall not tackle 
this problem, but rather limit ourselves to display these polynomials in the 
example at hand. They are 

A 2 - hiX +(hl-h 3 ) = 

/ 1 \ (6.21) 
H - h 2 X + [hih 2 - h 4 ) = 

The important idea emerging from the previous discussion is that the GZ equa- 
tions are often coupled with a special system of coordinates related with the 
Poisson pair. 

Exercise 6.2 Check that the polynomials ( |6.21|) define a system of Darboux- 
Nijenhuis coordinates for the pair (u , N) considered above. 



6.4 Separation of Variables 

We start from the classical Stackel theorem on the separability, in orthogo- 
nal coordinates, of the Hamilton-Jacobi equation associated with the natural 
Hamiltonian 

H(q, P ) = \Y^9 1 \<1)PI + Vfa, ...,q n ) (6.22) 

on the cotangent bundle of the configuration space. According to Stackel, 
this Hamiltonian is separable if and only if there exists as invertible matrix 
S(qx, ... ,q n ) and a vector U(qi, ... ,q n ) such that H is among the solutions 
(Hi, . . . , H n ) of the linear system 

1 n 

j=i 

and S and U verify the Stackel condition: 

The rows of S and U depend only on the corresponding coordinate. 

This means for instance that the elements Sij and U\ depend only on the first 
coordinate q±, and so on. Such a matrix S is called a Stackel matrix (and U a 
Stackel vector). 
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The strategy we shall follow to prove the separability of the Hamilton- Jacobi 
equations associated with the GZ vector fields X\ and X 3 on the manifold S4 
considered above, is to show that the Darboux-Nijenhuis coordinates allow to 
define a Stackel matrix for the corresponding Hamiltonians. 

The construction of the Stackel matrix starts from the matrix F which relates 
the vector field Xi and X 3 to the Nijenhuis tensor iV (see equation ( |6.15[ )). One 
can prove that this matrix satisfies the remarkable identity 

N*dF = FdF . (6.24) 

This is a matrix equation which must be interpreted as follows: dF is a matrix 
of 1-forms, and N* acts separately on each entry of this matrix; FdF denotes 
the matrix multiplication of the matrices F and dF, which amounts to linearly 



combine the 1-forms appearing in dF. In our example, equation ( 6.24 ) becomes 

N*d(TrN) = -d(detiV) + (TriV) d(TvN) 
N*d(detN) = + (detiV) d(TrJV) 



Exercise 6.3 Check that this equations are verified by the Nijenhuis tensor (|6.13|) 



We leave for a moment the particular case we are dealing with, and we suppose 
that, on a symplectic bi-Hamiltonian manifold fulfilling the conditions of Sub- 
section |6~3| , a family of n vector fields (Xl, X 3 , . . . , X 2n -i) is given. We assume 
that they are Hamiltonian with respect to Pq, say, X21-1 — PodCi, and that 
there exists a matrix F such that 

n 

NX 2i -i = F i X 2i-i for all i . (6.26) 



Finally, we suppose that F satisfies condition (|6.24j ). Then, from the matrix F 
we build up the matrix T whose rows are the left-eigenvectors of F. In other 
words, we construct a matrix T such that 

F = T _1 AT , (6.27) 

where A = diag(Ai, . . . , A n ) is the diagonal matrix of the eigenvalues of F, coin- 
ciding with the eigenvalues of N. The matrix T is normalized by imposing that 
in each row there is a constant component. A suitable normalization criterion, 
for instance, is to set the entries in the last column equal to 1. 

Theorem 6.4 // the matrix F verifies condition ( |6'. 2J$ (as it is always true in 
our class of examples), then the matrix T is a (generalized) Stackel matrix in 
the Darboux-Nijenhuis coordinates. 

This theorem means that the rows of the matrix T verify the following gen- 
eralized Stackel condition: The entries of the first row of T depend only on the 
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canonical pair (Ai,//i), those of the second row on (A 2 ,/z 2 ), and so on. With 
respect to the classical case recalled at the beginning of this lecture, we notice 
that by generalizing the class of Hamiltonians considered, we have been obliged 
to extend a little bit the notion of Stackel matrix. However, this extension 
does not affect the theorem of separability. Indeed, as a consequence of the 
fact that the matrix F is defined by the vector fields X 3 , . . . , X 2n _i) them- 
selves through equation ( 6.26|) , one can prove that T is a Stackel matrix for the 



corresponding Hamiltonians (C 1; . . . , C n ). 

Theorem 6.5 The column vector 

U = TC , (6.28) 

where C is the column vector of the Hamiltonians (Ci, . . . ,C n ), verifies the 
(generalized) Stackel condition in the Darboux-Nijenhuis coordinates. This means 
that the first component o/U depends only on the pair (Ai,/xi), the second on 
(A2, H2), and so on. 

We shall not prove these two theorems here, preferring to see them "at work" 
in the example at hand. First we consider the matrix T. Due to the form ( |6.15| ) 
of the matrix F, it is easily proved that 



T 



V A 2 I j 



(6.29) 



Indeed, the equation TF = AT follows directly from the characteristic equation 
for the tensor N. It should be noted that the matrix T has been computed 
without computing explicitly the eigenvalues Ai and A 2 . It is enough to use 
the first of equations (|6.21|) , defining the Darboux-Nijenhuis coordinates. The 
matrix T clearly possess the Stackel property (even in the classical, restricted 
sense) . 

The vector U can be computed as well without computing explicitly the 
coordinates (Xj,fij). It is sufficient, once again, to use the equations (|6.21|) . We 
now pass to prove that equation ( |6.28j ), in our example, has the particular form 

1 2 1 2 2 (6 ' 30) 
-/z 2 — -A 2 — E\ 2 = A 2 Ci + C 2 . 

We notice that proving this statement is tantamount to proving that the fol- 
lowing equality between polynomials, 

/i(A) 2 - A 5 = 2C(A) , (6.31) 

is verified in correspondence of the eigenvalues of N. This can be done as follows. 
Let us write the polynomials defining the Darboux-Nijenhuis coordinates in the 
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symbolic form 



A 2 = eiA + e 2 

(6.32) 

V = /iA + f 2 ■ 

The coefficients (ej, fj) of these polynomials must be regarded as known func- 
tions of the coordinates on the manifold. By squaring the second polynomial 
and by eliminating A 2 by means of the first equation, we get 

fi 2 = / 2 ( ei A + e 2 ) + 2A/ 2 A + hh = {flei + 2/i/ 2 )A + (f 2 e 2 + f 2 ) . (6.33) 

In the same way we obtain 



A 5 = A • A 4 = A[(e? + 2 ei e 2 )A + (e 2 e 2 + e 2 )] 

= (e 4 + 3e 2 e 2 + e 2 2 )X + [e\e 2 + 2e x e 2 2 ) 



(6.34) 



Finally, 

C(A) = C A 2 + d\ + C 2 = (C oei + d)A + (C e 2 + C 2 ) . (6.35) 



By inserting these expressions into equation Q6.31 ), we see that the resulting 



equation splits into two parts, according to the "surviving" powers of A: 

A : (f 2 ei + 2fJ 2 ) - (e 4 + 3e 2 e 2 + e 2 ) = 2(C Qei + C x ) 
1 : (e 2 + e 2 + e 2 ) - (e?e 2 + 2eie 2 ) = 2(C e 2 + C 2 ) . 

This method allows to reduce the proof of the separability of the Hamilton- 
Jacobi equation(s) to the procedure of checking that explicitly known functions 
identically coincide on the manifold. 

We end our discussion of the separability at this point. Our aim was simply 
to introduce the method of Poisson pairs, and to show by means of concrete 
examples how it can be profitably used to define and solve special classes of 
integrable Hamiltonian equations. We hope that the examples discussed in 
these lectures might be successful in giving at least a feeling of the nature and 
the potentialities of this method. 
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